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Abstract 

We consider within QCD collinear factorization the process p + p — > jet + jet + X, 
where two forward high-p^ jets are produced with a large separation in rapidity Ay 
(Mueller-Navelet jets). In this case the (calculable) hard part of the reaction receives 
large higher-order corrections ~ a" (Ay) n , which can be accounted for in the BFKL ap- 
proach with next-to-leading logarithmic accuracy, including contributions ~ a™ (Ay)™ . 
We calculate several observables related with this process, using the next-to-leading or- 
der jet vertices, recently calculated in the approximation of small aperture of the jet 
cone in the pseudorapidity-azimuthal angle plane. 
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1 Introduction 



We consider the inclusive production at high energies of two forward high-fc-p jets in proton- 
proton collisions, 

p(Pi) +P(p 2 ) -Met^jJ +jet(fc j2 ) +X , (1) 

which are detected in the fragmentation regions of two colliding protons, p(pi) and p(p2), 
and are separated by a large interval of rapidity Ay, the so-called Mueller-Navelet (MN) 
process pQ . This is considered as an important process for the manifestation of the BFKL [2] 
dynamics at hadron colliders, such as Tevatron and LHC. 

The theoretical description of this process is based on the QCD collinear factorization. 
Neglecting higher-twist contributions (terms suppressed with respect to the leading scaling 
asymptotic by additional inverse powers of the hard scale), the process can be viewed as 
started by two hadrons each emitting one parton, according to its parton distribution function 
(PDF), with the subsequent partonic hard scattering, see Fig.[TJ Collinear factorization allows 
to systematically resum the logarithms of the hard scale, calculating the standard DGLAP 
evolution [3] of the PDFs and the fixed-order radiative corrections to the parton scattering 
cross section. 

On the other side, in our kinematics at large squared center-of-mass energy y/s, when 
the rapidity gap between the two produced forward jets is large, the BFKL resummation of 
energy logarithms comes into play, since large logarithms of the energy compensate the small 
QCD coupling and must be resummed to all orders of perturbation theory. 

In comparison to the fixed-order DGLAP calculation, where an almost back-to-back emis- 
sion is expected, the BFKL calculation assumes more emission of partons between the two 
jets and leads generically to a larger cross-section and to a reduced azimuthal correlation 
between the detected two forward jets. 

At present the BFKL approach provides a general framework for the resummation of 
energy logarithms in the leading logarithmic approximation (LLA), which means resummation 
of all terms (a s ln(s)) n , and in the next-to-leading logarithmic approximation (NLA), which 
means resummation of all terms a s (a s ln(s)) n . Such resummation is process-independent and 
is encoded in the Green's function for the interaction of two Reggeized gluons. The Green's 
function is determined through the BFKL equation, which is an iterative integral equation, 
whose kernel is known at the next-to-leading order (NLO) both for forward scattering (i.e. 
for t = and color singlet in the t-channel) [U [5] and for any fixed (not growing with energy) 
momentum transfer t and any possible two-gluon color state in the t-channel 0, IE] • 

The process-dependent part of the information needed for constructing the cross section 
for the production of Mueller-Navelet jets is contained in the impact factors for the transition 
from the colliding parton to the forward jet (the so called "jet vertex"). Mueller-Navelet jet 
vertices were calculated with NLO accuracy in [9]. The results of [9] were then used in [TQl [TT| 
for a numerical estimation in the NLA of the cross section for Mueller-Navelet jets at LHC and 
for the analysis of the azimuthal correlation of the produced jets. This numerical analysis 
followed previous ones [121 [32] based on the inclusion of NLO effects only in the Green's 
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Figure 1: Mueller-Navelet jet production process. 



functions. Recently we performed a new calculation [2] of the jet impact factor, confirming 
the results of [9]. 

Although NLO jet vertices were obtained for a general jet algorithm, the implementation 
of these results into the cross section and the other Mueller-Navelet jet observables calculation 
requires a rather complicated numerical evaluation, see [10]. Recently the NLO impact factor 
for the production of forward jets was calculated by two of us [TJ)] in the "small-cone" approx- 
imation (SCA) [151 d7], i.e. for small jet cone aperture in the rapidity-azimuthal angle plane. 
The use of the SCA allowed to get a simple analytic result for the jet vertices in the so-called 
(u, n) -representation (the jet vertices projected on the eigenfunctions of the BFKL kernel) 
which can be easily implemented in numerical calculations. It is the aim of the present paper 
to obtain predictions for the Mueller-Navelet jet process cross section and for the azimuthal 
angle decorrelation observables in the SCA using the results of [T5] . 

The paper is organized as follows: in the next Section we discuss the kinematics and recall 
the formulae for the Mueller-Navelet jet process cross section; in Section 3 we present our 
results; Section 4 contains the discussion of results and our conclusions. 



2 Mueller-Navelet jet cross section 

It is convenient to define the Sudakov decomposition for the momenta of the jets; one has 



x Jx s 



Pi + k Jl± , 




Ji ' 



(2) 





(3) 



kj 2 



Xj 2 P2 + 



Pi + k J 2 1 , 



Xj 2 S 
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Here as power-suppressed correction we neglect the proton mass, p\ = p\ = 0, therefore p\ 
and p 2 are taken as Sudakov vectors satisfying 2{pip 2 ) = s. 

For the forward jets the longitudinal fractions x j x 2 ~ are related to the jet rapidities 
in the center-of-mass system by 

1, x \ s , dx J± 

2 k \ x * 

1 x\s dx j2 
y 2 = --\n-—, dy 2 = 

2 k% x j2 

so that the rapidity gap between the two jets is given by 



Ay = Y = In Z JlXj f . (4) 

I ^ Jl I I ^2 I 

For the Mueller-Navelet jet process, the rapidity gap Y is to be taken much larger than unity, 
thus implying the kinematics 

a > k\ 2 > A 2 QCD , (5) 
and the transverse momenta of the jets are assumed to be of similar order of magnitude, 
K Ji K .h- 

In QCD collinear factorization the cross section of the process ([T]) reads 

11 

dxj 1 dxj 2 d 2 kj 1 d 2 kj 2 . J 1 J 2 ^ 1 ^ F ^j( 2 ''^' F ^ dxj x dxj 2 d 2 kj 1 d 2 kj 2 ' ^ ^ 

where the i, j indices specify the parton types (quarks q = u,d,s,c,b; antiquarks q = 
u, d, s, c, b; or gluon g), fi(x, ftp) denotes the initial proton PDFs, the longitudinal fractions of 
the partons involved in the hard subprocess are Xx j2 , as shown in Fig. [U \ip is the factorization 
scale, d&ij(xix 2 s, hf) is the partonic cross section for the production of jets and s = x\x 2 s is 
the squared center-of-mass energy of the parton-parton collision subprocess. We use the MS 
scheme for the ultraviolet and collinear factorizations. 

At lowest order each jet is generated by a single parton having high transverse momentum 
and the partonic subprocess is given by an elementary two-to-two scattering. In the discussed 
Mueller-Navelet kinematics the higher-order contributions to the partonic cross section have 
to be resummed using BFKL approach. In the BFKL approach [2], the cross section of the 
hard subprocess reads 

d&i j(xiX 2 s) 1 f d 2 qi -> f d 2 q 2 



Vi^so.xukj^x^) / -^Vj{-q 2 ,SQ,x 2 \kj 2 ,xj 2 ) (7) 



dxj x dxj 2 d 2 kj Y d 2 kj 2 (2n) 2 J q 2 ' ' ' 1 J q 2 ' 

5+ioo 

f du ( xtx 2 s\ u 

x J 2^(^rJ G »^- 
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This representation for the cross section is valid with NLA accuracy. Here Vi(q\, s , x\\ kj x , XjA 
and Vj(—q 2 ,s ,X2]kj 2 ,xj 2 ) are the jet vertices (impact factors) describing the transitions 
parton i (xipi) — > jet (fcjj and parton j (x 2 p 2 ) — ^ jet (fcj 2 ), in the scattering off a Reggeized 
gluon with transverse momentum q\ and q 2 , respectively. The artificial scale so is introduced 
in the BFKL approach to perform the Mellin transform from the s-space to the complex 
angular momentum plane and cancels in the full expression for the cross section with the 
NLA accuracy. The Green's function in ([7]) obeys the BFKL equation 



<5 (<fi - = w G u {qx,q 2 ) - J d q K(q\,q) G u (q,q 2 ) , (8) 

where K(q\, q 2 ) is the BFKL kernel. 

In what follows we proceed along the lines similar to ones used in Ref. [18J. It is convenient 
to work in the transverse momentum representation, defined by 

Q \Qi) = QM} , (9) 

(9i|?2> = 5 (2) (<fi - q 2 ) , (A\B) = (A\k)(k\B) = j d 2 kA(k)B(k) ; (10) 

the kernel of the operator K is 

K(q 2 ,q l ) = {q 2 \K\q 1 ) (11) 
and the equation for the Green's function reads 

l = (u- K)G W , (12) 

its solution being 

G U] = {u-K)- 1 . (13) 

The kernel is given as an expansion in the strong coupling, 

K = a s k° + a\K x , (14) 

where 

a s = (15) 

7T 

and N c is the number of colors. In Eq. (fl4"]l K° is the BFKL kernel in the LLA, K 1 represents 
the NLA correction. 

To determine the partonic cross section with NLA accuracy we need an approximate 
solution of Eq. f lT3"j) . With the required accuracy this solution is 



G w = (w - c^T 1 + (u - ask )- 1 {a 2 s K 1 ^ (to - a s K°)- 1 + O (a^ 1 ' 

The basis of eigenfunctions of the LLA kernel, 
K°\n, v) = x(n, u)\n, v) , x(n, v) = 2-0(1) --0 77 + 7T + iv J ~ V> 77 + _ iu ) » ( 17 ) 



(16) 



2 2 / r \2 2 



is given by the following set of functions: 



(q\n,u) 



7T 



V2 



2\ w ~2 e in4 



18) 



here is the azimuthal angle of the vector q counted from some fixed direction in the transverse 
space, cos0 = q x /\q\- Then, the orthonormality condition takes the form 



(n', v'\n, v) 



d 2( l (->2\ i "- i v'- 1 i(n- 

2tt 2 {q > 



5{u — u') 5 n 



(19) 



The action of the full NLA BFKL kernel on these functions may be expressed as follows: 



K\n, v) = a s (n R )x(n, v) \n, v) + a 2 s (ii R ) [ x (1) (n, v) + T^x(n, v) ln(/4) ) \n, v) 



Po 



0o 



AN, 







dv 



(20) 



where \xr is the renormalization scale of the QCD coupling, the first term represents the 
action of LLA kernel, while the second and the third ones stand for the diagonal and the 
non-diagonal parts of the NLA kernel and we have used 



Po 



where rif is the number of active quark flavors. 



llNr 2n 



f 



(21) 



The function x i n > v )i calculated in [19] (see also |20j), is conveniently represented in the 
form 



0o 



X 2 (n, v) - ^x(n, v) - ix'{n, v) ) + x(n, v) 



(22) 



where 



+ 



7r 2 sinh(7ri/) 
2 v cosh 2 ('7ri/ > ) 



7T 2 -4 



X {n, v) - 6C(3) - X "{n, v) + 2 <f>(n, v) + 2 <j>(n, -v) (23) 



3+1 + -^ 



n f \ 11 + 12z/ 2 



N% J 16(1 + v 2 ) 



>n0 



1+^ 



1 + 4t/ 2 



0(n, v) 



dx 



x 



-l/2+iu+n/2 



1 + X 



1 V +1 



iV c 3 J 32(1 + v 2 ) 
C(2) ) +Li 2 (x) + Li 2 (-x 



»2 



(24) 



E 

fc=0 



+ lnx [ ^(n + 1) -^(1) +ln(l + x) + 
l)k +i 



-x) 



OO i. 



k=l 



k + n 



E 



^ (A; + n) 5 



'1 



fc + (n + l)/2 + iV 



[i>'(k + n + 1) - ^>'(fc + 1) + (-l) fc+1 (/3'(A; + n + 1) + /?'(*; + 1)) 



k + (n + l)/2 + iv 



11/ 



, fz + 1 



-11/ 



Li 2 (x) 



dt 



ln(l - t) 



Here and below x'( n i v ) = dx( n , is)/dis and x"( n i v ) = d 2 x( n i is)/d 2 is. 

For the quark and the gluon jet vertices in (J2J) the projection onto the eigenfunctions of 
LO BFKL kernel, i.e. the transfer to the (v, n)-representation, is done as follows: 



V(q 



i +00 

r~ = S / dis^i(is,n)(n,u\q 1 ) 



V(-q2) 



+00 
+00 „ 

^2 / dv$ 2 (v,n)(q 2 \n,v) 



$i(i/,n) = / d 2 q l 



,2,. ^(<fl) 1 ^-2\^-| ^in^x 



d%^^^{q 2 2 r u ^ 2 ^ 2 



The vertices can be represented as an expansion in a s , 

^1,2(1, f) = a s (fi R )v ly2 (n, v) + a^(/i fi )^(n, 1/) . 



(25) 



(26) 



In Eqs. (125]) and ( 126]) we suppressed for brevity the partonic indices i,j and the other argu- 
ments in t> 1)2 . The explicit forms of LLA and NLA jet vertices in the (is, n)-representation 
both for the quark and gluon cases can be found in [15] . In particular for the LLA quark 
vertices one has 



vi(n,i/) 



v 2 \n,v) 



hl^\r- 3/2 ^s( XJl - Xl ) 



2 ./^(^)-^-3/2 e -m(^ 2+7 r) 5( 



x Jl ~ x 2j 



(27) 



where Ca = N c , Cf = (N 2 — 1)/2N C , the angle 0j 2 + n in the last equation appears due to 
the fact that the Reggeon momentum which enters the second vertex is — q 2 . Note that in 
LLA vertices the partonic and the jet longitudinal momentum fractions coincide. 

The partonic cross section can be written with NLA accuracy as follows 

_L~> +°° 



da(x\x 2 s) 



dxj x dxj^d 2 kj x d 2 kj 2 (2ir) 



X\X 2 S 

so 



xaldiRjv^n, is)v 2 (n, is) 



«s(Mfl)x(™,^) 



+a s (n R )\n 



X\X 2 S 



( v[ \n,is) v { 2 \n,is) 

1 + a s (n R ) — r- + — 

y v t (n, is) v 2 (n, is) 

0o , . r ,10 



X(n,u) + ^X(n, v) 



-X(n, is) + 



+i 



dl n ( VlM 

\v 2 (n,v) 

dis 



+ 21n/4 



(28) 
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For the subsequent calculation it is convenient to make the substitution 



X\X 2 S 

so 



■* J Ji x ,hS \ / x l \ (-^2 



«0 



£j 2 



(29) 



and to assign the last two factors in the r.h.s. to the corresponding jet vertices. This procedure 
affects only the NLA parts of the jet vertices, since for the LLA vertices Xj = xj t . Also with 
NLA accuracy, one can make in (1281) the replacement 



1„ _> l n 



SO 



(30) 



This procedure allows to perform in the MN-jet cross section first the integration over partonic 
momentum fractions, before taking the sum over n and the integration over v\ it allows also 
to consider together the contributions of quarks and gluons to the jet vertices. 



The differential cross section has the form 

da 1 



C + 2 cos(n0) C n 



71=1 



(31) 



where 



>Ji - VJ2 

r-2lT 1-2-K 



— 7r, and 



C m = d$ Jx \ d(f)j 2 cos[m(0 Jl - <pj 2 - tt)] 
Jo Jo 



da 



dy.hdy j2 dlfcjjd^jjd^ja 



(32) 



In particular, taking into account the Jacobian of the transformation from the variables 
kji, to the variables \kj t \, yj t , and the ^-dependence of LLA jet vertices, see (j2*7| . we get 



r 



\kj x \\k 



+00 



dv 



J 2 | J -°° 



% Ji% J2S 
so 



X 



Xa2(/ifl)c 1 (n, v, |A-j 1 |,x Jl )c 2 (n, v, \kj 2 \,x j2 ) 
, .(^\n,v, IfejJ.sjJ 4 1} (n,z/, |fcj 2 |,x j2 ) 

1 + a s {fI R ) I 13 1 =; 

ci(n,i/, I^jJ^jJ 02(71, z/, |/cj 2 |,x j2 ) 



-2, , 1 / x Ji x J 2 s \ I -1 \ , A> / \ / / \ . 10 , 1 A*b 



+a 2 (/i R )ln 



s 



Vi'Va. 



where 



c 1 (n,z/,|fc|,x) = 2,/^(fc 2 r- 1 / 2 [^/ B (x lA iF)+X;/ 9 (^A*F 



C 



c 2 {n,u, \k\,x) 



ci(n, v, \k\,x) 



a=q,q 



Ci(n,v, \k\,x) 



TT 



1 ic f 



F t2 



-1/2 fd( 
7 



c 



-as(M«)x(™^) 



(33) 



(34) 
(35) 
(36) 
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( P«(0 + ^P 9q (C)) In K ~ In R {P w (0 + P gq (0} ~ % ^ ^(1 - C) 

+C A S(1 - (x(n, 7 )ln|l + ^ + y + i (V (l + 7 + 5) " <l>' (\ - l) - x\n, 7 ))) 

+ (1 + e) \c A ( ^gg^ - c- (^) ) + ( Q 



+ 2 



ln(l - C) 

i-C 



7T 2 9\ 5n/ 



C, 21nC 

C 2 2 i-C 



5(l-C)(C F (31n2- y --l - g 



+C A ( + CVC + (ca-Ji + 2C A In ^ + CVC 27 (x(™, 7) - 2 In ()) 



c c 



(^(0 + 2n f ^P qg (C)) ln^ - 2r 2 ^ln J R(P ra (C) + 2n,P w (0) - ^ In -£1,5(1 - C) 
+ C A 5(1 - C) (x(n, 7 )ln^ + l + ^ + l(^(l + 7 +^)-^Q- 7 )- X 2 (n, 7))) 
+ 2^(l-C^)(Q-2 + C c)lnC>^f^) 



+ 71/ 



l + (r^o;- 2+cc ; 



(1 + C 27 )x(n, 7) -2 In C+^/2 



C_ 2 

C 2 ' 



2CC g + (C 2 + C 2 ) (gx(n, 7) + ^3) - ^(1 - C) 



4 1} (n,z/, |A;|,x) 



cS 1} (n, 1/, |fc|,x) 



(37) 



Here ( = 1 — (,'y = iu — 1/2, Pij(() are leading order DGLAP kernels. For the ^1,2,3 functions 
we have the results: 



j _c 2 



c 



^(1,1+7-22 + 7-U) 2^(1,1 + 7 + 1,2 + 7 + 1,0 



f-7-1 



+ 7 + 1 



(38) 



( _ 2 .. ^(1, -7 - 2 1 - 7 - 2 C) 2 Pi(l, -7 + f , 1 - 7 + f , 



i + 7 



2 ' 



+ (1 + C 27 ) (x(n, 7 )-21nC)+21nC] , 



lnC+^f^(x(n,7)-21nC) 



(39) 
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Figure 2: Y dependence of Co for 



h J2 I 



35 GeV at yfs = 14 TeV. 



C r C 



1HC+ 1 5 7 (x(n iT )-21nC) 



(40) 



The factor ^-M^kO^) appears in (13T)j) due to extra-contributions attributed to the jet ver- 
tices, as discussed after Eq. ( 129|) . Note that the C n coefficients do not depend on the azimuthal 
angles of the jets, <pj 1 and 0j 2 , 

(41) 



C n = C n [ yj^yjz, k Ju kj 2 ,n R , s 



they depend instead on the jet rapidities, the transverse momenta and on the factorization, 
renormalization and energy scale parameters. 



An alternative way to present the differential cross section, equivalent to the formula (|33 
in the NLA, is the so-called exponentiated form (see Ref. [IB]), 



C 



exp 



I k Ji 1 1 k Ji I 



dv exp 



(Y - Y ) (a s (fi R )x(n, v) 



(42) 
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/4 



h 2 h 2 



xaf(/iji)ci(n,i/, 1^1,2:^,^)02(71,1/, \kj 2 \,x j2 ,ii F ) 
a ^ ^ I cS 1} (n, t/, IfcjJ^j^/if) + 4 1} (n, t/, \kj 2 \,xj 2 ,fi F ) 

Qi(7l, I/, IfcjJjXju/iF) C 2 (n, Z/, IfcjjIjXjaj/XF) 
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Table 1: Values of Co in the LLA, in the NLA and in the NLA with LO impact factors for 
l^jj — \kj a \ = 35 GeV at y/s = 14 TeV, corresponding to the data points in Fig. |2j The 

optimal values of Y and ur — [irI \J l^jJI^jal f° r C^ hA ^ are given in the fourth and fifth 
columns, while those for (j^ LA / LO IF ) are given in the last two columns. 
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Table 2: Values of C in the LLA and in the NLA for l&jj = |jfej 2 | = 20 GeV at y/s = 14 TeV, 
corresponding to the data points in Fig. [3j The last two columns give the optimal values of 



Y and n R = /jl r / \J {kj^kj^. 



where Y is defined in Eq. (pE} and we have introduced the variable 



V n = ln 



so 



\ k Ji\\ k J2\ 



(43) 



Below we will discuss the differential cross section integrated over the jet azimuthal angles 



C = J dfadfada , 

the coefficients C n and the moments of the azimuthal decorrelations, which are defined as 

(cos(n^)) = f = — . 44 

J dcp^dcp^da C 



3 Results 



In this section we present our results for the dependence on Y of the functions C n . In what 
follows we take the factorization and renormalization scales equal to each other, //p = /jr. We 



10 



: 

100 r 



10 



0.1 



+ 



Vs = 14 TeV 



+ LLA 
x NLA 



J_ 



J_ 



k = k T = 20 GeV 



6 7 8 9 10 Y 

Figure 3: Y dependence of Co for \kj^ = \kj 2 \ =20 GeV at y/s = 14 TeV. 



perform our calculation both in the LLA and in the NLA. In the former case, the expression 
for C„ reads 



C 



LLA 



da LL A 



dyj^dyj^kj^ d\k 



J*\ 



\k 



h -h I 



dv exp 



(Y - Y )a s {n R )x(n, v) 



(45) 



xa s (fi R )ci(n, v, \k Jl \,x Jl )c 2 (n,v, \kj 2 \,xj 2 ) . 
For our NLA analysis we use the exponentiated representation given in Eq. fj42]h 

For the center-of-mass energy \fs we take the LHC design value 14 TeV. We fix the jet 
cone size at the value R = 0.5, in order to compare our predictions with the forthcoming LHC 
data. We study Mueller-Navelet jets with symmetric and asymmetric values of the transverse 
momenta, in particular, consider the choices: j/cjj = \kj 2 \ = 35 GeV, = \kj 2 \ = 20 GeV 
for the symmetric cases and jfcjj = 20 GeV, \kj 2 \ = 35 GeV for asymmetric option. 

Moreover, to make possible the comparison with the experiments at present LHC energy, 
we perform calculations for yfs = 7 TeV, where we consider the = \kj 2 \ = 35 GeV case. 



Following a quite recent CMS study [2T] , we restrict the rapidities of the Mueller-Navelet 
jets to the region 3 < \yj\ < 5. We will present our results for Co, i.e. the differential cross 
section integrated over the jet azimuthal angles, the coefficients C n for ti^O, and (cos(n0)) 
versus the relative rapidity, Y = yj 1 — yj 2 . For our choice of forward jets rapidities, Y 
takes values between 6 and 10. Our approach is similar to the one used in [10], we introduce 
rapidity bins with step Ayj equal to 0.5, so the considered values for jet rapidities and rapidity 
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Figure 4: Y dependence of Cq for j&jj = 20 GeV and |fcj 2 | = 35 GeV at yfs = 14 TeV. 
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Table 3: Values of C in the LLA and in the NLA for \k. h \ = 20 GeV and \k, h \ = 35 GeV 
at y/s = 14 TeV, corresponding to the data points in Fig. HI The last two columns give the 

optimal values of Y and n R = /j, R / \J\ k j x \ \ k j 2 \ . 
difference are 

{{y.h)i} = {3.0,3.5,4.0,4.5,5.0} 
{{VM = {-3.0,-3.5,-4.0,-4.5,-5.0} 

and {Yi} = {6.0, 6.5, 7.0, 7.5, 8.0, 8.5, 9.0, 9.5, 10.0}. Then we evaluate the following sum 

CniY^^VSCniiyj^^iyj^^-Yi) 
where the sum runs over the possible values of (yjjj for a given Y{. 

In our analysis we use the PDF set MSTW2008nnlo [22] and the two-loop running coupling 
with a s (M z ) = 0.11707. 
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Figure 5: Y dependence of Cq for \kj t \ = \kj 2 \ = 35 GeV at y/s = 7 TeV. 



Our predictions depend on the values of energy and renormalization scales, s and fi R . 
For the analysis in the LLA we fixed the values of these scales, and s , as suggested by 
the kinematics of the process, i.e. fi 2 R = sq = \ k j x 1 1 k j 2 \ . 

In general LLA results depend very strongly on sq and fin, and one really needs to proceed 
to the NLA analysis in order to reduce this scale dependence and to have some reliable 
predictions for observables. One should stress that the dependence of the correlations C n on 
the scales and s Q cancels with NLA accuracy; nevertheless in both representations ( 13"3|) 
and there unavoidably exist contributions subleading to NLA, depending on fi R and s Q , 
whose numerical impact is important for the considered kinematics, therefore we need some 
prescription for the choice of these scales. 
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Table 4: Values of Cq in the LLA and in the NLA for jfcjj = \kj 2 \ = 35 GeV at \fs = 7 TeV, 
corresponding to the data points in Fig. [5J The last two columns give the optimal values of 



Y and n R = UrN \kj x \\kj 2 \. 
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Figure 6: Y dependence of C\/Cq (left) and C2/C0 (right) for [kj^ = \kj 2 \ = 35 GeV at 
y/s = 14 TeV. Please note that LLA data at Y = 8 and 10 from the present work and from 
Ref. [10] are barely distinguishable in these plots. 
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Table 5: Values of Ci/Cq = (cos0) in the LLA, in the NLA and in the NLA with LO impact 
factors for = \kj 2 \ = 35 GeV at yfs = 14 TeV, corresponding to the data points in 

Fig. [Srieft). The optimal values of Y and n R = fj, R / \\kj 2 \ for c[ NLA ^ /Cq NLA ^ are given 
in the fourth and fifth columns, while those for (j( NLA / LO if)^q(nla/lo if) given in the 
last two columns. 



Following Ref. [IB], we use here an adaptation of the principle of minimal sensitivity 
(PMS) [23], which consists in taking as optimal choices for fi R and So those values for which 
the physical observable under examination exhibits the minimal sensitivity to changes of both 
of these scales. The motivation of this procedure is that the complete resummation of the 
perturbative series would not depend on the scales fi R and so, so the optimization method is 
supposed to mimic the effect of the most relevant unknown subleading terms. 

In our search for optimal values, we took integer values for Y in the range 0-5 and values 
for jj, R given as integer multiples of y |fcji||^j 2 |, 



HR = n R \j\kj 1 \\kj 2 \, (46) 

taking the integer n R in the range 1-15. The systematic uncertainty of the optimization 
procedure in the determination of observables, which will be discussed below, originates from 
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Table 6: Values of C2/C0 = (cos<p) in the LLA, in the NLA and in the NLA with LO impact 
factors for = \kj 2 \ = 35 GeV at y/s = 14 TeV, corresponding to the data points in 

Fig. E]4eft). The optimal values of Y and tir = /j,r/ \J\kj 1 \\kj 2 | for (^ nla )/q( nla ) are g[ Yen 
in the fourth and fifth columns, while those for ^ NLA / LO if)^^nla/lo if) given in the 
last two columns. 
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Figure 7: Y dependence of C\/Cq (left) and C2/C0 (right) for l&jj = \kj 2 \ = 20 GeV at 
= 14 TeV. 



the resolution of the grid in the ur - Yq plane. This uncertainty has been estimated as the 
standard deviation of the optimal value from the determinations in the nearest neighbors of 
the grid. The error bars around the NLA data points presented in the figures below represent 
this uncertainty. We did not evaluate the impact on our predictions of the PDF uncertainties, 
since we expect it to be of the same size as determined in Ref . [10] , where the same PDF set 
adopted here is used. 

Let us start with the cross section integrated over the jet azimuthal angles, Cq. We found 
that for this observable a stationary point in the - Y plane could always be singled out, 
typically a local maximum. For y/s = 14 TeV our results, in [ G "ya ] units, are presented in 
Figs. |2rll] and in Tables dH31 results for y/s = 7 TeV are given in Fig. and Table SJ The 



optimal values of Y and % = y l^/JI^I are also reported in the tables. As in previous 
works [IS] , the optimal values of the energy scales turn to be far from the kinematic scale. On 
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Table 7: Values of Cx/Cq = (cos0) in the LLA and in the NLA for \kj x \ = \kj 2 \ = 20 GeV 
at -\/s = 14 TeV, corresponding to the data points in Fig. [TJ^left). The last two columns give 

the optimal values of Y and tir = [Mr/ y [kj^k^ \. 
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Table 8: Values of C 2 /C = (cos(20)) in the LLA and in the NLA for \k^ \ = \kj 2 \ = 20 GeV 
at y/s = 14 TeV, corresponding to the data points in Fig. fright). The last two columns give 

the optimal values of Y and tir = /j,r/ yikjJlkj^. 

the other hand, the uncertainty related of our optimization procedure, described above, turns 
out to be small, therefore our NLA results for the cross section integrated over jet azimuthal 
angles, presented in Figs. [2H51 have relatively small "error bars" . 

The other issue we addressed in a similar manner is the analysis of the observables Ci/Cq 
and C2/C0, which encode the first two nontrivial angular decorrelations: (cos 0) and (cos(20)) . 
For \kj x \ = \kj 2 \ = 35 GeV at y/s = 14 TeV our results are presented in Fig. [6] and Tables |5] 
and El For the smaller values of jet transverse momenta, \kj^\ = \kj 2 \ = 20 GeV, they are 
given in Fig. [7] and Tables [7] and El Finally, for the asymmetric case, = 20 GeV and 
\kj 2 \ = 35 GeV, we present our results for (cos0) and (cos(20)) in Fig. [8] and Tables l9j fTDl 

In our search for optimal scale values we found that these observables, angular decorrela- 
tions, are less stable than Co, thus making our analysis more complicated. Indeed, only for a 
few values of Y; t we could find a stationary point. In other cases we found a local maximum 
only in the direction of one of the two parameters. When this happened, we took as "optimal" 
value of the observable of interest the one exhibiting the least standard deviation from values 
taken in the nearest neighboring points in the adopted grid. As a result the uncertainties 
related with our optimization procedure are larger for the angular decorrelations than for 
the integrated cross sections discussed above, therefore Figs. EHH1 show larger "error bars" in 
comparison to the NLA results presented in Figs. [2HI3 Note that for the asymmetric case, 
\kj^\ = 20 GeV and \kj 2 \ = 35 GeV, the search for optimal scale values was more problem- 
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Figure 8: Y dependence of C\/Cq (left) and C2/C0 (right) for 
GeV at v/i = 14 TeV. 



20 GeV and \k 



■h\ 



35 



Y 


£,(LLA) ^(LLA) 


£,(NLA) ^(NLA) 


Y 


n R 


6 


0.173 


0.555(13) 


1 


8 


7 


0.129 


0.530(25) 


1 


5 


8 


0.094 


0.475(17) 


2 


5 


9 


0.068 


0.387(11) 


1 


8 


10 


0.0481 


0.3477(95) 


1 


9 



Table 9: Values of Ci/C = (cos</>) in the LLA and in the NLA for \k^\ = 20 GeV and 
\kj 2 \ = 35 GeV at y/s = 14 TeV, corresponding to the data points in Fig. IHl^left). The last 



two columns give the optimal values of Y and ur = /xr/v |&jJ \kj 2 \. 



atic for the observable (cos(20)}, which is affected by uncertainties larger than those of the 
correlation (cos0). 

Even if in general the energy dependence of the cross section and the azimuthal decorre- 
lations is driven mainly by the kernel, in the considered kinematics the contribution of the 
NLO corrections to impact factors happened to be important. In order to show this, in the 
analysis at |/cjj = \kj 2 \ = 35 GeV and y/s = 14 TeV we calculated the coefficient C and the 
correlations C\/Cq and C2/C0 using the NLA BFKL kernel together with the LO impact fac- 
tors (properly modified with the inclusion of the NLO terms which guarantee the cancellation 
of the hr- and so-dependence in the NLA expressions for the coefficients C n , see Ref. |24j). 
We find that the NLO corrections to impact factors are relevant especially at large values of 
Y, as shown in Figs. H] and EJ 

In the following Section we discuss the results presented here, make some comparison with 
Refs. pm |TT] and draw our conclusions. 
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Table 10: Values of C 2 /C = (cos(20)) in the LLA and in the NLA for (fcjj = 20 GeV and 
\kj 2 \ = 35 GeV at y/s = 14 TeV, corresponding to the data points in Fig. fright). The last 

two columns give the optimal values of Y and tir = /ir/ \J\kjJ\kj^\. 

4 Discussion 



In this paper we considered in NLA BFKL approach the Mueller-Navelet jet production 
in proton-proton collisions, using the results for NLA jet vertices obtained recently in the 
"small-cone" approximation. Having a simple analytic result for the jet vertices, projected 
on the eigenfunction of LLA BFKL equation ((is, n)-representation), one can implement them 
easily in the calculation of the Mueller-Navelet jet cross section. All necessary formulae are 
presented in Section 2. 

We confirm the observation found earlier in the works devoted to forward electroproduc- 
tion of a pair of vector mesons [IB] , and to Mueller-Navelet jets production [ID], that NLA 
corrections to the impact factors (jet vertices) are very important and can not be ignored in 
a consistent NLA BFKL analysis. 

Our numerical results, presented in Section 3, depend on the energy and renormalization 
scales, so and hr. The dependence of the coefficients C n on these scales cancels with NLA 
accuracy after the inclusion of NLA corrections to jet vertices. Nevertheless, due to next-to- 
NLA contributions depending on hr and so, the observables we calculated are sensitive to 
the choice of these scales. Here, following Ref. [18J, we used an optimization procedure, based 
on the principle of minimal sensitivity [23], which consists in taking as optimal choices for 
fiR and So those values for which the physical observable exhibits the minimal sensitivity to 
changes of both these scales. 

The small-cone approximation, which we adopted here, is expected to be an adequate 
tool. Indeed, it is known that in the general case the dependence of the cross section on the 
jet cone parameter has, in the limit R — > 0, the form da ~ A In R + B + 0(R 2 ) (see, for 
instance, [T6] and Appendix C there). Indeed, in SCA the coefficients A and B are evaluated 
exactly. The neglected pieces 0(R 2 ) for typical R values are presumably less important than 
the other uncertainties of our NLA BFKL calculation, in particular those related with the 
choice of the scales fiR and so, which mimic in our method the effect of the most relevant 
unknown next-to-NLA BFKL terms. 

To support this statement it seems natural to make a comparison of our results obtained 
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in SCA with the numerics presented in [TO] ITT] , where the jet cone size was treated exactly. 
We would make such comparison for \kj x \ = \kj 2 \ = 35 GeV and yfs = 7 and 14 TeV. For this 
purpose we present in Table [TTJ our NLA results obtained with the above-discussed optimal 
scales setting (presented in the third, fifth and seventh columns) and compare them with 
those for Co, (cos0) and (cos(20)) taken from Tables 1, 5 and 9 of Ref. [.10] and reported in 
the second, fourth and sixth columns of our Table [TT]) . Moreover, in Table [121 we show our 
NLA results in the case when kinematic values of the scales were used, fj, R — s — Ikj^kj^. 

Let us discuss the numbers presented in Tables [TTJ and O Firstly one needs to say that 
NLA results obtained with our formulae at kinematic scale setting, fi R — so — |&JiH&j 2 | ; 
can not be regarded as acceptable predictions. Indeed we got in this case results for (cos</>) 
(third column of Table [T2|) which are away from the kinematic range of the cosine function 
and for the integrated cross section Co (second column of Table [T2"|) we obtained even a 
negative value in the case of the largest rapidity difference, Y = 10. This is related to the fact 
that NLO corrections to the jet vertices are negative and very large in absolute value when 
the kinematic scale setting is used. A similar observation was done in Refs. [T8] . where the 
electroproduction of a pair of vector mesons was considered. In such a situation we actually 
need the PMS procedure in order to make reliable predictions. 

Returning to the comparison with the results of JTU] , we observe that our results when the 
kinematic scales were used (second, third and fourth column of Table IT2"|) are quite different 
from those of Ref. [10] (second, fourth and sixth column of Table: ITT]) , where the same 
scales were used. As discussed above, in our opinion such big difference cannot be attributed 
only to the SCA. Most probably the real source of discrepancy is in the actually different 
representations of the NLA BFKL cross sections adopted in the two calculations, which are 
equivalent with NLA accuracy. This issue is not simple and definitely deserves a further study. 
At this stage we can only compare our results obtained with the PMS procedure with those 
of HO]. 

For the integrated cross section, Coflthe results presented in the second and third columns 
of Table [TT] are rather close each other. Let us discuss now the observable (cos 0) . Our 
predictions for it in the case of = \kj 2 \ = 35 GeV and a/s = 14 TeV are shown in Fig. O 
In Table [TT] we compare them (fourth and fifth columns) with those obtained in Ref. |10] . 
In this case we see more difference between two approaches, in particular our results show a 
clear tendency for (cos0) to decrease with Y, whereas pjj] predicts a flat V-dependence of 
(cos0). Instead, for the observable (cos(20)) the agreement between two approaches turns to 
be rather good, as sees in Table [TTJ sixth and seventh columns. 

Similarly, we compare our predictions for Co at the present LHC energy a/s = 7 TeV with 
those of Ref. [TT]. The agreement is fair, except for the case Y = 10, as shown in Table [T3T 

Traditionally the BFKL predictions are assumed to be compared with the fixed order 

1 Our LLA results coincide with those of [W\ with high accuracy. We note, in passing, that the results 
quoted in Tables 15, 19, 35, 38, 53 and 56 of Ref. [TP] , giving the coefficients C\ and C2 for several values of 
the jet kinematics, should be multiplied by a factor two to correctly reproduce the values of the ratios C1/C0 
and C2/C0 quoted in other tables of that paper. We stress that this normalization problem does not affect 
any of the comparisons presented in this work between our results and those of Ref. 10 . 



19 



Y 




^y(NLA) 


\cos <p) 


^COS (p) 


{COS{Z(p ) ) 


{cos{Z(p ) ) 




fRpf nrni 


Vi prp 
lid t. 


[Rpf nrni 

[rcei. jiujj 


11"! "I 




Vi pro 
11"! " 


6 


0.606 


0.613 


0.851 


0.7124 


0.512 


0.457 


7 


0.670 


0.675 










8 


0.289 


0.299 


0.777 


0.608 


0.383 


0.338 


9 


0.0474 


0.0496 










10 


0.00238 


0.00257 


0.753 


0.525 


0.317 


0.290 



Table 11: Comparison of our predictions for the observables Co, (cos0) and (cos(20)) in the 
NLA for \k 3l \ = \k Ja \ = 35 GeV at y/s = U TeV with those of Ref. [10]. 
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Table 12: Values of Co, (cos0) and (cos(20)) in the NLA for jfcjj = \kj 2 \ = 35 GeV at 
y/s = 14 TeV for the kinematic values of the scales, fj, R = Sq = {kj^kj^ or Y = and 
n R = 1. 

DGLAP ones, trying to find a kinematic range where possible experiment can discriminate 
between these two approaches. For Mueller-Navelet process the relevant parameter which 
can describe the separation of these two regimes is 77 = a s (\kj\)Y, which has the meaning 
of the mean number of hard undetected partons inclusively produced in the process. For 
\kj\ = 35 GeV and Y = 6 4- 10 this parameter takes the values rj = 0.82 ~ 1.37. In Section 3 
we presented also our predictions for smaller jets transverse momenta, \kj\ =20 GeV, where 
7] = 0.92 4 1.53 for Y = 6 -j- 10. In this case the BFKL description is expected to give results 
more different with respect to the NLO DGLAP ones. We hope that experiments with such 
Mueller-Navelet jet transverse momenta will be possible in the future at LHC. 
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Table 13: Comparison of our predictions for the observable Co for \kj t \ = \kj 2 \ = 35 GeV at 
y/s = 7 TeV with those of Ref. HQ. 
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Our PMS optimization procedure for the kinematics considered gives us optimal energy 
and factorization scales values which are substantially larger than the scale given by the 
kinematics, [i 2 R = \ k j 1 \ \ k j 2 \ . This fact indicates the presence of important contributions 
subleading to the NLA. Therefore the estimates of the uncertainties in our predictions should 
be taken with care. Nevertheless they reflect the reliability of PMS method in the considered 
kinematics. Note that despite very large negative contributions to the MN-jet cross section 
coming both from the NLA corrections to the BFKL kernel and NLA corrections to the jet 
vertices with respect to LLA MN-jet cross section, we got with our PMS procedure rather 
precise results in all cases except for the (cos(20)) observable at [fcjj = 20 GeV and \kj 2 \ = 
35 GeV and y'i = 14 TeV. It would be very interesting to confront our predictions with 
experiment. 
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